A Single Trapped Ion as a Time-Dependent Harmonic Oscillator 
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We show how a single trapped ion may be used to test a variety of important physical mod- 
els realized as time-dependent harmonic oscillators. The ion itself functions as its own motional 
detector through laser-induced electronic transitions. Alsing et al. [Phys. Rev. Lett. 94, 220401 
(2005)] proposed that an exponentially decaying trap frequency could be used to simulate (thermal) 
Gibbons-Hawking radiation in an expanding universe, but the Hamiltonian used was incorrect. We 
apply our general solution to this experimental proposal, correcting the result for a single ion and 
showing that while the actual spectrum is different from the Gibbons-Hawking case, it nevertheless 
shares an important experimental signature with this result. 
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I. INTRODUCTION 

The time-dependent quantum harmonic oscillator has 
long served as a paradigm for nonadiabatic time- 
dependent Hamiltonian systems and has been applied to 
a wide range of physical problems by choosing the mass, 
the frequency, or both, to be time-dependent. The ear- 
liest application is to squeezed state generation in quan- 
tum optics P, 0, [I[ , in which the effect of a second-order 
optical nonlinearity on a single-mode field can be mod- 
eled by a harmonic oscillator with a frequency that is 
harmonically modulated at twice the bare oscillator fre- 
quency. It was subsequently shown that any modulation 
of the frequency could produce squeezing [J] , and thus the 
same model could be used to approximately describe the 
generation of photons in a cavity with a time-dependent 
boundary [1, Q. 

The model has been used in a number of quantum 
cosmological models. In Ref. Q, a time-dependent fre- 
quency has been used to explain entropy production in a 
quantum mini-superspace model. The model, with both 
mass and frequency time-dependent, has been particu- 
larly important in developing an understanding of how 
quantum fluctuations in a scalar field can drive classical 
metric fluctuation during inflation [1, Q. In a cosmo- 
logical setting the time-dependence is not harmonic and 
is usually exponential. In all physical applications, of 
course, the model is only an approximation to the true 
physics, and its validity can be tested only with consid- 
erable difficulty, especially in the cosmological setting. 
Here we propose a realistic experimental context in which 
the time-dependent quantum harmonic oscillator can be 
studied directly. 

Many decades of effort to refine spectroscopic measure- 
ments for time standards now enable a single ion to be 
confined in three dimensions, its vibrational motion re- 
stricted effectively to one dimension, and the ion cooled 



to the vibrational ground state with a probability greater 
than 99% flol |. Laser cooling is based on the ability to 
couple an internal electronic transition to the vibrational 
motion of the ion [ll|. These methods can easily be ex- 
tended to more than one ion and their collective normal 
modes of vibration [l^] • Indeed so carefully can the cou- 
pling between the electronic and vibrational states be 
engineered that is is possible to realise simple quantum 
information processing tasks [H, [3| • We use the control 
of trapping potential afforded by ion traps, together with 
the ability to reach quantum limited motion, to propose 
a simple experimental test of quantum harmonic oscilla- 
tors with time-dependent frequencies. We also make use 
of the ability to make highly efficient qua ntum measure- 
ments, based on fluorescent shelving [10j, to propose a 
practical means to test our predictions. 

In this paper, we calculate the excitation probability 
of a trapped ion in a general time-dependent potential. 
When beginning in the vibrational ground state of the 
unchirped trap and starting the chirping process adia- 
batically, the excitation probability is simply related to 
the Fourier transform of the solution of the Heisenberg 
equations of motion (which is also the same as the trajec- 
tory of the equivalent classical oscillator). We compare 
our result with that of Ref. [l5| for the case of a single 
ion undergoing an exponential frequency chirp. The cited 
work attempts to use this experimental setup to model 
a massless scalar field during an inflating (i.e., de Sitter) 
universe, which would give a thermal excitation spectrum 
as a function of the detector response frequency [1 61 ] . The 
analysis is incorrect, however, because the wrong Hamil- 
tonian was used. Nevertheless, the corrected calculation 
presented here also gives an excitation spectrum with 
a thermal signature, although the particular functional 
form is different. 



II. GENERAL SOLUTION 
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The quantum Hamiltonian for a single ion in a time- 
dependent harmonic trap can be well-approximated in 
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one dimension by 
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where v(t) is time-dependent but always assumed to 
be much slower than the timescale of the micromo- 
tion For emphasis, we have indicated the explicit 
time-dependence of the frequency z/; we will often omit 
this from now on. Working in the Heisenberg picture, we 
get the following equations of motion for q and p: 



P 

q= M' 

p = -Mv 2 q . 



(2) 
(3) 



Dots indicate total derivatives with respect to time. Dif- 
ferentiating again and plugging in these results gives 



= q + v 2 q , 

i> , 2 
= p — 2—p + v p 
v 



(4) 
(5) 



As we shall see, only Eq. is necessary for calculating 
excitation probabilities, so we will focus only on it. These 
equations are operator equations, but they are identical 
to the classical equations of motion for the analogous 
classical system. Interpreting them as such, we will la- 
bel the two linearly independent c-number solutions as 
h(t) and g(t), where the following initial conditions are 
satisfied: 



M0) = <?(0) = 1 



and 



h(0) = .g(0) = , (6) 



Writing q(0) = qo and p(0) = po, the unique solution for 
q to the initial value problem above is 



q(t) = q h{t) 



Po 
M 



9(t) 



(7) 



By differentiating and using the relations above, we know 
also that 



p{t) = Mq h(t) + p g(t) 



(8) 



To check our math, we can verify that [q(t),p(t)] = ih, 
which is fulfilled if and only if the Wronskian W(h, g) of 
the two solutions is one for all times — specifically, 



hg — hg = I 



(9) 



where we have assumed that [go,Po] = 

Moreover, if the initial state at t — is symmetric with 
respect to phase-space rotations, then we have additional 
rotational freedom in choosing the initial quadratures. 
(This would be the case, for instance, if we start in the 
instantaneous ground state.) Notice that Eq. ([7]) can be 
written as the inner product of two vectors: 
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qo, 
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h(t),v g(t) 



(10) 



(and similarly for Eq. |5])), where we have normalized 
the quadrature operators to have the same units. As an 
inner product, this expression is invariant under simulta- 
neous rotations of both vectors. Thus, if the initial state 
possesses rotational symmetry in the phase plane, then 
the rotated quadratures are equally as valid as the orig- 
inal ones for representing the initial state, which means 
that an arbitrary rotation can be applied to the second 
vector above without changing any measurable property 
of the system. This freedom can be used, for instance, to 
define new functions h'(t) and g'(t) that are more con- 
venient for calculations, where the linear transformation 
between them and the original ones (with prefactors as 
in Eq. (jTDJ) ) is a rotation. We will use this freedom in the 
next section. 

One reason why ion traps have become a leading im- 
plementation for quantum information processing is the 
ability to efficiently read out the internal electronic state 
using a fluorescence shelving scheme • As the internal 
state can become correlated with the vibrational motion 
of the ion, this scheme can be configured as a way to 
measure the vibrational state directly [l7} . To correlate 
the internal electronic state with the motion of the ion, 
an external laser can be used to drive an electronic tran- 
sition between two levels \g) and |e), separated in energy 
by Hloa- The interaction between an external classical 
laser field and the ion is described, in the dipole and 
rotating-wave approximation, by the interaction-picture 
Hamiltonian Oil 



Hl = — ifif2( 



i6q(t) 
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i8q(t) 



(11) 



where fio is the Rabi frequency for the laser-atom inter- 
action, u>l is the laser frequency, k is the magnitude of 
the wave vector k, which makes an angle 9 with the trap 
axis, q(t) is given in Eq. ([7]), and 



a± (t) — e a± 



(12) 



The electronic-state raising and lowering operators are 
defined as a + = \e)(g\ and er_ = |g)(e|, respectively, and 



A = UJA - UJL 



(13) 



is the detuning of the laser below the atomic transi- 
tion. We can construct a meaningful quantity that char- 
acterizes the "size" of q(t) based on the width of the 
ground-state wave packet for an oscillator with frequency 
v(t), namely yj h/2Mv(t). As long as this quantity is 
much smaller than k cos throughout the chirping pro- 
cess, then we can expand the exponentials in Eq. (111[) to 
first order and define the interaction Hamiltonian Hi be- 
tween the electronic states and vibrational motion (still 
in the interaction picture) by 



Hi = hflok cos 8q(t) (e 



-iAt. 



+iAt 



(14) 



where we have assumed that u>l is far off-resonance, and 
thus A ^ 0. 
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Using first-order time-dependent perturbation theory, 
the probability to find the ion in the excited state is 



p{1) = JpJ o dt iJ o dt 2 (H I (t 1 )V e H I (t 2 )) 



flgk COS 



T pT 

dt\ I dt 2 e 

n JO 



lA(tl - t2) (q(h)q(t 2 )) , 

(15) 



where V e = l V ib ® l e )(e| is the projector onto the ex- 
cited electronic state (and the identity on the vibrational 
subspace). We always assume that the ion begins in the 
electronic ground state. If the ion also starts out in the 
instantaneous vibrational ground state for a static trap 
of frequency v$ — v(0) at t — (which is most useful 
when the chirping begins in the adiabatic regime), then 
we can evaluate the two-time correlation function as 



<9(ti)g(*a)> groU nd = (lo) Kh)h(t 2 



M 2 
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[h{t x )g{t 2 )-h{t 2 )g{t x )] 
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h 
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h(ti) - iv g{t\) h(t 2 ) + iv g{t 2 ) 
f(ti)f*(t 2 ) , 



(16) 



where we have used the facts that for the vibrational 
ground state, (q^) = ((po/Mv ) 2 ) = hj2Mv a and 

(qoPo) = | {{qo,Pq} + [9o,Po]) = and we nav e de- 

fined the complex function 



f(t) = h(t) - iu g(t) 



(17) 



which is the solution to Eq. (0| with initial the conditions, 
/(0) = 1 and /(0) = -iv . Plugging this into Eq. (HSJ) 
gives, quite simply, 



P« - (« r?o) 2 |F| 2 , 



where 



F = 



dtt 



iAt 



/(*) 



(18) 



(19) 



and we have defined the unitless, time-dependent Lamb- 
Dicke parameter [l(j as 



n(t) 



Ihk 2 cos 2 6 
2Mv{t) 



(20) 



and T70 = ??(0). Recalling that f(t) can be considered a 
complex c-number solution to the equations of motion for 
the equivalent classical Hamiltonian, Eq. (fTS)) shows that 
the excitation probability is simply related to the Fourier 
transform of the classical trajectories when beginning in 
the vibrational ground state. 



III. EXPONENTIAL CHIRPING 

Recent work [l5} has suggested that an exponen- 
tially decaying trap frequency has the same effect on 
the phonon modes of a string of ions as an expand- 
ing (i.e., de Sitter) spacetime does on a one-dimensional 
scalar field [HI]. An inertial detector that responds to 
such an expanding scalar field would register a thermal 
bath of particles, called Gibbons-Hawking radiation 
Ref. [l5| suggests that the acoustic analog [l9[ of this 
radiation could be seen in an ion trap, causing each ion 
to be excited with a thermal spectrum with temperature 
hK,/2TrkB, as a function of the detuning A, where k is 
the trap-frequency decay rate. The analysis used an in- 
correct Hamiltonian that neglected squeezing and source 
terms that have no analog in the expanding scalar field 
model but which are present when considering trapped 
ions in this way, and the results are incorrect. In this sec- 
tion, we revisit this problem and calculate the excitation 
probability for a single ion in an exponentially decaying 
harmonic potential, as a function of the detu ning A. 

We write the time-dependent frequency as [20l | 



v{t) = v e 



This results in 



q + vle- 2Kt q = 



Solutions with initial conditions © are 



hit) = ^ 
V ' 2k 



(21) 
(22) 

(23) 
(24) 



where the time dependence is carried in v — v{t) from 
Eq. (|2"Tj) , and J n and Y n are Bessel functions. We could 
plug these directly into the formulas from the last section, 
but we will simplify the calculations by considering the 
limits of slow and long-time frequency decay, represented 
by 



V » K 



and 



u e kT <C k , 



(25) 



respectively. This allows us to do several things. First, it 
allows us to use the usual ground state of the unchirped 
trap at frequency isq as a good approximation to the 
ground state of the expanding trap at t — 0, since at that 
time the system is being chirped adiabatically. This is 
important because it allows the experiment to begin with 
a static potential, which is useful for cooling. Second, it 
allows us to simplify h(t) and g(t) using the phase-space 
rotation freedom discussed above. Using asymptotic ap- 
proximations for the Bessel functions in the coefficients, 



J ' (?) * 



2k 
ttv 



■ cos 



K 



(26) 
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we get 
h{t) 



2k 



sin^lo (-J + cosi/3 Jo y^j , (28) 



7T^o 

2k 



cos <p Yc 



sin </3 Jq 



(29) 



where </j = vq/k — ■n/A. Since we are taking the initial 
state to be the ground state, which is symmetric with 
respect to phase-space rotations, we can use the freedom 
discussed in the previous section to undo the rotation 
represented by Eqs. (|2"5|) and (|29[) and define the simpler 
functions 



i(<W(t) 



(30) 
(31) 



The primes are unnecessary due to the symmetry of the 
initial state, so we drop them from now on and plug di- 
rectly into Eq. (JTTJ) : 



*(;)-<) 



(32) 



where H„ is a Hankel function of the first kind. The 
integral in Eq. lfH?|) can be evaluated in the limits (j2"5|) 
using techniques similar to those used in Ref. [15J. First, 
define 

e a = — . t = a — nt , it = e T , and x — A/k . 

K 

(33) 

The integral in question then becomes (neglecting the 
prefactor) 



dt e 



iAt tt( 



- lAt H^\e a " Kt ) 



>'C) ■ 

i / dTe- M(a - r) ^ 1) (e T ) 

K Ja-nT 
-ixa poo 

" dr e lXT H { Q l \e T ) 



K 

-ixa poo 



(34) 



Inserting a convergence factor with x — > x — ie, and then 
taking the limit e — > + , we can use the formula 



duu^H^iu) 



T{ix/2) 



(e™ - l)T(l-ix/2) 



(35) 



to evaluate 



, ,9 ttmi 1 
F = — - — 

M 2k k 2 



T(ix/2) 



2T7VQ 



r(i -ix/2) 
i 



K 3 X 2 (gTrx _ ]_)2 



(36) 



When plugging in for the dummy variables (I33p . this 
gives 



pW = (fl %) : 



, 27Wn 



1 



kA 2 (e 



tA/k 



I) 2 



(37) 



The calculated result from Ref. [15j for a single ion is 



P 



(i) 



GH 



(fio»7o) 



2tt 



1 



kA e 27rA / K - 1 



(38) 



which contains a Planck factor with Gibbons- 
Hawking [16j | temperature T — hK/2itkB but is 
different from the actual result for a single ion, given by 

Eq. dnzD 

Several things should be noted about these functions. 
First, they both break down as A — ► because of the ap- 
proximation made in obtaining Eq. (|14p . They also fail 
if the time-dependent Lamb-Dicke parameter ([20]) ever 
becomes too large throughout the chirping process. Fur- 
thermore, most cases of interest will be A ~ v§ (the first 
red sideband) and near A ~ — 1/ (the first blue side- 
band), which means that |A| ^> k, since u$ 3> k. The 
first red sideband represents a detector that requires the 
absorption of one phonon (plus one laser photon) in order 
to excite the atom — the usual thing we mean by "particle 
detector" when the particles are phonons. The first blue 
sideband, on the other hand, represents a detector that 
emits a phonon in order to excite the atom (along with 
absorbing one laser photon). 

There are a couple of ways to compare these functions. 
First, we can take the ratio of the two for both the red- 
and blue-sideband cases. In both cases, we obtain 



P 



(i) 



GH 



^L(l + 2e-^) 



(39) 



plus terms of order 0(e~ 27r ' A '/ K ). Since |A| ~ v , the 
prefactor is close to one, and the second term is very 
small (since 3> k). Furthermore, it is cumbersome 
to directly compare the measured probability to the full 
function (with all the prefactors). It is often easier in- 
stead to make measurements on both the first red side- 
band and the first blue sideband and then take the ratio 
of the two. The constant prefactors disappear in this 
calculation, and both functions then have the same ex- 
perimental signature: 



P«(A) _ P«(A) 



PW(-A) P«(-A) 



_ -27tA/k 



(40) 



which is that of a thermal distribution with tempera- 
ture T = h,K/2itkB, which is of the Gibbons-Hawking 
form [If} with the expansion rate given by k. There- 
fore, although the Hamiltonian used in the calculations 
in Ref. [l5| was missing terms, the intuition (at least for 
a single ion) was correct in that the actual experimental 
signature in this case matches that of an ion undergoing 
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thermal motion in a static trap, where the temperature 
is proportional to k. 

To see whether this experiment is feasible, we must ex- 
amine the validity of our approximations. For a typical 
trap, we expect that isq ~ 1 MHz, and thus if we take 
K ~ 1 kHZ, we easily satisfy the first of conditions (|25l) . 
namely i>q 3> k. The second of these conditions gives a 
constraint on the modulation time T. For these param- 
eters we expect that T ~ a few msec. This is compat- 
ible with typical cooling and readout time scales and is 
less than those for heating due to fluctuating patch po- 
tentials [nj. Thus, this is a realizable experiment with 
current technology. 

IV. CONCLUSION 

We have shown that a single trapped ion in a modu- 
lated trapping potential can serve as an experimentally 
accessible implementation of a quantum harmonic oscilla- 
tor with time-dependent frequency, including robust con- 
trol over state preparation, manipulation, and measure- 



ment. The ion itself serves both as the oscillating particle 
and as the local detector of vibrational motion via cou- 
pling to internal electronic states by an external laser. 
For the case of a general time-dependent trap frequency, 
we calculated the first-order excitation probability for the 
ion in terms of the solution to the classical equations of 
motion for the equivalent classical oscillator. We applied 
this general result to the case of exponential chirping and 
corrected the calculation in Ref. [l5| for a single ion. We 
found that while the results from the two calculations dif- 
fer, the experimental signature in both cases is the same 
and equivalent to that of a thermal ion in a static trap. 
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